Reduced Thermal Conductivity of Nanowires and Nanoribbons with Dynamically Rough Surfaces 
and the "Problem of One-Dimensional Heat Conductors" 
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We present analytical model and molecular dynamics simulations of phonon heat transport in nanowires 
and nanoribbons with anharmonic lattices and dynamically rough surfaces and edges. In agreement with re- 
cent experiments on heat transport in single-crystalline silicon nanowires with rough surfaces, our model and 
simulations predict finite and length-independent phonon thermal conductivity in such quasi-one-dimensional 
systems, in contrast to anomalous phonon thermal conductivity of corresponding momentum-conserving sys- 
tems with atomically smooth surfaces, divergent with the system length. Within our model, the main cause 
of thermal conductivity reduction is momentum-nonconserving scattering of longitudinal acoustic phonons by 
anharmonic side phonon leads in quasi-one-dimensional phonon waveguide with dynamically rough surface or 
edge layers. 
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Thermal conductivity (TC) of low-dimensional nanostruc- 
tures has recently attracted much interest in connection with 
the search for technologically feasible materials for thermo- 
electric applications [1, 2]. High-efficiency thermoelectric 
materials are important for power-generation devices, which 
are designed to convert waste heat into electrical energy, and 
therefore they are expected to play an increasingly important 
role in meeting the energy challenge of the future. The re- 
duction of the TC without significant reduction of the elec- 
tric conductivity is the most important strategy for enhanc- 
ing the thermoelectric figure of merit, the parameter which 
measures the material potential for thermoelectric applica- 
tions [1, 2]. Recent experiments have shown that TC of sili- 
con nanowires (NWs) can be dramatically reduced by surface 
roughness [3, 4]. Furthermore, the phonon contribution to TC 
of a single-crystalline Si NW with diameter D ^ 50 nm was 
shown to approach the limit of amorphous Si and therefore 
a single-crystalline Si NW with such small diameter behaves 
like a "phonon glass" [3]. This observation cannot be ex- 
plained by the existing theories, although the important role 
of surface roughness in phonon scattering in NWs has been 
emphasized [5, 6, 7, 8]. Molecular dynamics (MD) modeling 
of diamond nanorods with functionalized surfaces also shows 
significant TC reduction in such low-dimensional systems [9]. 

On the other hand, quasi-ballistic long wavelength acous- 
tic phonons in disordered harmonic or ideal anharmonic one- 
dimensional (ID) systems give rise to TC coefficient, diver- 
gent with the system length, see, e.g., [10, 11, 12, 13, 14, 15, 
16]. As it was shown for the first time by Peierls, only the 
momentum-nonconserving scattering of phonons can limit TC 
of the system [17]. In this Letter we present analytical model 
and MD simulations of the TC of NWs and nanoribbons 
(NRs) with dynamically rough surfaces and edges which show 
dramatic decrease of TC in such partially-disordered systems, 
in contrast to TC of corresponding ideal systems (with atom- 
ically smooth surfaces and edges). The analytical model de- 
scribes ID lattice of coupled oscillators in which each oscil- 
lator is weakly coupled to weakly-anharmonic side oscilla- 
tory chain, or side phonon lead (SPL), which models dynami- 



cal surface roughness of the quasi-one-dimensional phonon 
waveguide (QIDPW), see Ref. [18]. Within this model, 
the main cause of TC reduction is momentum-nonconserving 
scattering of longitudinal acoustic phonons by SPLs in dy- 
namically rough surface or edge layers of the QIDPW. Strong 
scattering of coherent longitudinal acoustic phonons in the 
QIDPW is a result of energy losses caused by the excitation 
of acoustic phonons in SPLs, which propagate only outwards 
and do not return back to the QIDPW in the form of coher- 
ent waves [18]. This occurs when the length of a SPL exceeds 
acoustic phonon mean free path (localization length) / 
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to anharmonic scattering in the lead (which gives I 
The effective length of the atomic or molecular SPL increases 
by random intersections with other SPLs in the rough surface 
or edge layer. 

The proposed phonon scattering mechanism gives finite and 
frequency-independent mean free path Iph of the propagat- 
ing along QIDPW longitudinal acoustic phonons in a wide 
frequency range, including low frequencies [18]. The value 
of phonon mean free path due to the dynamical-roughness- 
induced momentum-nonconserving scattering is given by the 
lattice spacing in the QIDPW a and coupling parameter be- 
tween a SPL and QIDPW g: Iph^^a/g. Coupling parame- 
ter g is determined by the ratio of the cross-section area of 
the rough layer [r^nDS) with the average root-mean-square 
roughness height S and cross-section area of NW with diame- 
ter D: gr^S/D [18]. For atomically-thick rough surface layer 
with (5>a, acoustic phonon mean free path lph^D{a/S) can 
become shorter than NW diameter [18]. Therefore this chan- 
nel of surface phonon scattering substantially reduces mean 
free path of "long longitudinal waves" which give dominant 
contribution to phonon TC of quasi-lD systems, see Refs. 
[10, 11, 16]. This in turn causes a substantial reduction of 
phonon TC in NWs and NRs with dynamically rough sur- 
faces or edges, which explains the significant TC reduction, 
down to the level of TC of amorphous Si, observed in single- 
crystalline Si NWs with rough surfaces [3]. The predicted 
reduction factor a/S is quantitatively consistent with the ob- 



2 



served reduction of TC of Si NWs with rough surface layers 
with controlled thicknesses S>a [3] with respect to TC of Si 
NWs with the same diameters and smooth surfaces (with i5~a) 
[18]. 

This mechanism of momentum-nonconserving acoustic 
phonon scattering gives finite and length-independent coeffi- 
cient of TC of QlDPWs with atomically rough surfaces or 
edges with 5>a, in contrast to the divergent with the sys- 
tem length coefficient of TC of anharmonic QlDPWs with 
atomically smooth surfaces or edges with (5<Ca, when the 
momentum-nonconserving channel of acoustic phonon scat- 
tering is closed and mean free path of long-wave phonons 
substantially increases, Ip/^ocuj^'^ , and becomes Ipii^D. The 
considered mechanism of acoustic phonon scattering by vi- 
bration of atomic or molecular chains in dynamically rough 
surfaces and edges of NWs and NRs is qualitatively differ- 
ent from usually considered acoustic wave scattering by static 
(geometric) roughness on stress-free surfaces of solids [5, 14]. 

Analytical model of phonon transport in quasi-lD 
phonon waveguides with side phonon leads. We consider 
ID lattice with period a of coupled oscillators, with coupling 
constant c and mass m, in which each oscillator is weakly cou- 
pled to SPL ("lateral" atomic chain with weakly-anharmonic 
coupling), with (linear) interatomic coupling ci=gc^c and 
atomic mass 'mi=gm-^m. Distribution of such SPLs along 
the QIDPW models its dynamical surface roughness [18]. 
The most important assumption of this model is that acous- 
tic phonon, propagating along the QIDPW, excites acoustic 
waves in SPLs which propagate along the leads only outwards 
the QIDPW and do not return back to the waveguide in the 
form of coherent waves. With this assumption, we obtain dis- 
persion equation for the dimensionless complex wave number 
fc||a of the damped acoustic phonons, propagating along the 
QIDPW [18]: 
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This equation predicts that phonons with frequencies lo 3> 
Lo* = g^Jc/m propagate quasi-ballistically, when 



k\\a 



(2) 



while low-frequency phonons, with u^oj*, propagate dif- 
fusively, when Lo—~iDphk'^\, with the diffusion coefficient 

Dph=a^ \/ cj mj g=Vphlph, where Yph=a^ elm and lpj^=a/g 
are velocity and energy mean free path of long-wave acous- 
tic phonons, cf. Eq. (2). Both quasi-ballistic and diffusive 
acoustic phonons contribute to phonon TC coefficient K,ph 
of QIDPW as Kph^CphDph^CphVpha/ g, where Cph is the 
contribution of longitudinal phonons to the specific heat of 
NW or NR. This predicts, in agreement with the observations 
[3, 4] and our MD simulations, the significant reduction of 
phonon mean free path and TC caused by phonon scatter- 
ing by dynamical surface roughness in NWs and NRs (with 
gr^5/D and (5>a). 
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Figure 1 : (a) Model of 2D atomic ribbon with rough edges which 
consists of M=\2 chains with ideal core of A/o=4 chains and two 
rough edges of M\=A chains withp=0.3 randomly deleted atoms, (b) 
Model of 3D nanowire with rough surface layers of 2 atom thickness, 
(c) Ideal 3D core. 



MD simulations of TC of 2D nanoribbons. We consider 
the atomic structure consisting of AI parallel chains of atoms 
placed in one plane. We consider the Hamiltonian of the scalar 
model of the ribbon lattice in which only the longitudinal 
displacements of atoms are taken into account and displace- 
ment of the (m, n)th atom from its equilibrium position. The 
Hamiltonian, which accounts for interaction between only the 
nearest-neighbor atoms (with unit mass): 
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where N is number of atoms in each chain, Vmn and Umn 
are potentials of the intra- and interchain interaction between 
(m, n) and {m, n + 1) atoms, and between {m, n) and (m + 
1, n) atoms, respectively. For the purpose of calculating the 
heat flux along the NR, we can define the corresponding total 
energy of the nth cross-section e„ and local energy flux j„ in 
the ribbon, which satisfy the continuity condition e„=7„-i — 
j„, see, e.g., Ref. [19]. 

We consider a NR built of M = Mq + 2Mi chains. To 
model two dynamically rough edges with widths Mi, we ran- 
domly delete some atoms in the chains m ~ 1 , . . . , Mi and 
TO = Ml + Mo + 1, Ml + Mo + Ml. Let < p < 1 be 
the probability of atom removal. In result of the random atom 
removal from the edge layers, some atoms in the edges will 
become completely isolated and should be deleted as well. 
Model NR with two rough edges with p ~ 0.3 randomly 
deleted atoms is shown in Fig. 1(a). 

Assuming for the certainty that U{p) = ^V{p), where 
V{p) = V,nn{p) and U{p) = Umn{p), wc consider V{p) = 
pV2, V{p) = p^/2 + pV4, V{p) = exp(-p) + p - 1, 
V{p) ~ [cxp(— p) — l]^/2 and V{p) = 1 — cos(p) as, respec- 
tively, the harmonic, Fermi-Past-Ulam (FPU), Toda, Morse 
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and rotational potentials. It is worth mentioning that the rib- 
bons with the harmonic, FPU, Toda or Morse interatomic po- 
tentials and ideal atomically smooth edges have infinite coef- 
ficients of TC in the limit of ^ cx). In purely harmonic sys- 
tems, acoustic phonons do not interact and there is no energy 
scattering in phonon thermal transport. Infinite TC coeffi- 
cients in ribbons with the anharmonic FPU, Toda or Morse in- 
teratomic potentials are related with the quasi-ballistic trans- 
port of the long-wave acoustic phonons with long mean free 
paths Ipfi oc Ideal ribbons with the rotational potential 

have finite TC in the limit of — > oo and finite temperatures 
[20,21]. 

The ribbon TC is found in two independent ways. The 
first method relies on direct modeling of heat transport along 
the ribbon. For this purpose, we consider a ribbon with to- 
tal length N + 2Na (Nq = 40 is the length of ribbon ends), 
which ends are placed in Langevin thermostats with different 
temperatures of the left, T+ = I.IT, and right,r_ = 0.9T, 
ends, where T is average temperature). In the middle part 
of the ribbon with length N the stationary heat flux J„ = J 
and linear temperature gradient of ribbon temperature T„ are 
established (where J„ is determined by the time-average of 
the local heat flux (jVi) in the stationary conditions). Thus the 
TC coefficient k{N) of the finite-length ribbon can be reliably 
determined as 



k{N) 
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where Mq is the width of central ideal strip, 2AIiq is total 
width of rough edges with filling fraction q=l — p < 1. TC 
coefficient k is determined from k{N) as the limit of ^ 
oo. The ribbon temperature profile T„ depends on a particular 
realization of the edge roughness and therefore it is necessary 
to perform averaging over its independent realizations. In our 
simulations, we averaged over 120 independent realizations 
of the edge roughness. 

We modeled heat transport in rough-edge NRs with cen- 
tral part length of A^=20,40,80, 160,320,640. Dependence of 
TC coefficients on the length of NR with rough edges, with 
p = 0.3 randomly deleted atoms, is shown in Fig. 2(a) for 
different interatomic potentials. As one can see in this figure, 
the rough-edge ribbon has finite coefficient of TC in the limit 
of infinite length for all the considered nonlinear interatomic 
potentials. In contrast to that, TC coefficient k{N) of the 
rough-edge ribbon with purely harmonic lattice continuosly 
decreases for N-^oo, see Fig. 2(b), and the limiting value of 
TC coefficient of such NR is zero, k(oo)=0. Figure 2(b) also 
shows that for the given width of the ideal ribbon core Mq, TC 
coefficient of the NR decreases with the increase of the rough 
edges width M i . Therefore the infinite-length ribbon or wire 
with purely harmonic lattice and dynamically rough edges can 
be considered as an ideal thermal insulator. 

We relate the reason for such strong difference between 
TC coefficients of rough-edge NRs with harmonic and an- 
harmonic interatomic potentials with the properties of SPLs 
in corresponding systems. In the anharmonic systems, long 
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Figure 2: (a) Dependence of TC coefficient k on length N of 2D 
ribbons with smooth and rough edges for the Toda (curves 1 and 4, 
temperature T = 1), Morse (curves 2 and 5, temperature T = 0.1) 
and FPU (curves 3 and 6, T=10) interatomic potentials, (b) Depen- 
dence of TC coefficient k on length A'^ of 2D ribbons with purely 
harmonic lattices and rough edge widths A/i=l and i\/i=4 (curves 
7 and 8) with p=0.3 randomly deleted atoms. Ideal core width is 
A/o=2, temperature T=l. 



acoustic waves, excited in SPLs, do not return back and the 
effective "internal radiative losses" result in finite phonon 
mean free path and phonon TC in the QIDPW, see Eqs. 
(1) and (2). In purely harmonic systems, the long acoustic 
waves, excited in SPLs, return back into the QIDPW and 
strongly suppress the transmission of low-frequency acous- 
tic phonons through the QIDPW because of destructive in- 
terference with the phonons propagating along the waveguide 
core (the phonon Fano resonance, see Ref. [18] for a brief 
review). This produces the effective stop band (or Anderson- 
Fano localization) for low-frequency acoustic phonons in the 
QIDPW with harmonic lattice and dynamically rough edges 
[18], and, correspondingly, the decrease of TC coefficient of 
the NR in the limit of N-^oo. Indeed, we observe only the 
localized vibration eigenstates in harmonic phonon spectrum 
of such system (not shown). 

In the approach based on Green-Kubo method, TC is de- 
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Figure 3: Dependence of TC coefficient of nanowires witli ideal 
(curve 1) and rough (curves 2 and 3) surfaces on length A'^, for the 
Toda (curves 1 and 2) and harmonic (curve 3) potentials. 

fined as integral of the autocorrelation of heat fluxes: 

t 



where J{t)=J2n J" (^) '■^^ average heat flux in the NR. Val- 
ues of k{N) obtained with the use of two different approaches, 
Eqs.(4) and (5), coincide with good accuracy. 

We also perform MD simulations of the temperature de- 
pendence of TC of rough-edge NRs with the FPU, Toda and 
rotational interatomic potentials (not shown). All the obtained 
dependencies confirm the finite coefficients of TC of the NRs 
in the low temperature limit T ^ when the coefficients of 
classical TC of corresponding NRs with atomically smooth 
edges diverge. 

MD simulations of TC of 3D nanowires. We consider 3D 
model system made of M x K parallel molecular chains as it is 
shown in Figs. 1(b) and 1(c). The Hamiltonian of the system 
we take as a 3D generalization of the Hamiltonian (3) of 2D 
ribbon scalar model. In Fig. 3 we show the length dependence 
of nanowire TC coefficient for the Toda and harmonic inter- 
atomic potentials. It demonstrates that similar to the case of 
2D nanoribbon, the dynamical roughness of the surface layer 
changes the increasing (divergent) with the length TC of 3D 
nanowire to the finite or decreasing with the nanowire length 
TC of 3D nanowire with, respectively, the anharmonic or har- 
monic interatomic potential. 

In summary, we present analytical model and molecular 
dynamics simulations of phonon heat transport in nanowires 



and nanoribbons with anharmonic lattices and dynamically 
rough surfaces and edges. In agreement with recent experi- 
ments on heat transport in single-crystalline silicon nanowires 
with rough surfaces, our model and simulations predict fi- 
nite and length-independent phonon thermal conductivity in 
such quasi-one-dimensional systems, in contrast to anomalous 
phonon thermal conductivity of corresponding momentum- 
conserving systems with atomically smooth surfaces, diver- 
gent with the system length. We also present the thermal- 
insulator-like heat transport in long-length nanowires and 
nanoribbons with purely harmonic lattices and dynamically 
rough surfaces, caused by the Anderson-Fano phonon local- 
ization. 
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